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Abstract. In the 2-dimensional curved 3-body problem, we prove 
the existence of Lagrangian and Eulerian homographic orbits, and 
provide their complete classification in the case of equal masses. 
We also show that the only non-homothetic hyperbolic Eulerian 
solutions are the hyperbolic Eulerian relative equilibria, a result 
that proves their instability. 



1. Introduction 

We consider the 3-body problem in spaces of constant curvature {k, ^ 
0), which we will call the curved 3-body problem, to distinguish it from 
its classical Euclidean (k = 0) analogue. The study of this problem 
might help us understand the nature of the physical space. Gauss 
allegedly tried to determine the nature of space by measuring the angles 
of a triangle formed by the peaks of three mountains. Even if the goal 
of his topographic measurements was different from what anecdotical 
history attributes to him (see [6]), this method of deciding the nature 
of space remains valid for astronomical distances. But since we cannot 
measure the angles of cosmic triangles, we could alternatively check 
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whether specific (potentially observable) motions of celestial bodies 
occur in spaces of negative, zero, or positive curvature, respectively. 

In [2], we showed that while Lagrangian orbits (rotating equilateral 
triangles having the bodies at their vertices) of non-equal masses are 
known to occur for k = 0, they must have equal masses for k ^ 0. 
Since Lagrangian solutions of non-equal masses exist in our solar system 
(for example, the triangle formed by the Sun, Jupiter, and the Trojan 
asteroids), we can conclude that, if assumed to have constant curvature, 
the physical space is Euclidean for distances of the order 10^ AU. The 
discovery of new orbits of the curved 3-body problem, as defined here in 
the spirit of an old tradition, might help us extend our understanding 
of space to larger scales. 

This tradition started in the 1830s, when Bolyai and Lobachevsky 
proposed a curved 2-hody problem, which was broadly studied (see most 
of the 77 references in P]). But until recently nobody extended the 
problem beyond two bodies. The newest results occur in P], a paper 
in which we obtained a unified framework that offers the equations of 
motion of the curved n-body problem for any n > 2 and k ^ 0. We also 
proved the existence of several classes of relative equilibria, including 
the Lagrangian orbits mentioned above. Relative equilibria are orbits 
for which the configuration of the system remains congruent with itself 
for all time, i.e. the distances between any two bodies are constant 
during the motion. 

So far, the only other existing paper on the curved n-body problem, 
treated in a unified context, deals with singularities, P], a subject we 
will not approach here. But relative equilibria can be put in a broader 
perspective. They are also the object of Saari's conjecture (see [7], |1]), 
which we partially solved for the curved n-body problem, [2]. Saari's 
conjecture has recently generated a lot of interest in classical celestial 
mechanics (see the references in [3], [5]) and is still unsolved for n > 3. 
Moreover, it led to the formulation of Saari's homographic conjecture, 
[7], |5], a problem that is directly related to the purpose of this research. 

We study here certain solutions that are more general than relative 
equilibria, namely orbits for which the configuration of the system re- 
mains similar with itself. In this class of solutions, the relative distances 
between particles may change proportionally during the motion, i.e. the 
size of the system could vary, though its shape remains the same. We 
will call these solutions homographic, in agreement with the classical 
terminology, [8]. 

In the classical Newtonian case, [8] , as well as in more general classi- 
cal contexts, P], the standard concept for understanding homographic 
solutions is that of central configuration. This notion, however, seems 
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to have no meaningful analogue in spaces of constant curvature, there- 
fore we had to come up with a new approach. 

Unlike in Euclidean space, homographic orbits are not planar, unless 
they are relative equilibria. In the case k > 0, for instance, the inter- 
section between a plane and a sphere is a circle, but the configuration 
of a solution confined to a circle cannot expand or contract and remain 
similar to itself. Therefore the study of homographic solutions that 
are not relative equilibria is apparently more complicated than in the 
classical case, in which all homographic orbits are planar. 

We focus here on three types of homographic solutions. The first, 
which we call Lagrangian, form an equilateral triangle at every time 
instant. We ask that the plane of this triangle be always orthogonal 
to the rotation axis. This assumption seems to be natural because, 
as proved in |2j, Lagrangian relative equilibria, which are particular 
homographic Lagrangian orbits, obey this property. We prove the ex- 
istence of homographic Lagrangian orbits in Section 3, and provide 
their complete classification in the case of equal masses in Section 4, 
for K > 0, and Section 5, for k < 0. Moreover, we show in Section 6 
that Lagrangian solutions with non-equal masses don't exist. 

We then study another type of homographic solutions of the curved 
3-body problem, which we call Eulerian, in analogy with the classical 
case that refers to bodies confined to a rotating straight line. At every 
time instant, the bodies of an Eulerian homographic orbit are on a 
(possibly) rotating geodesic. In Section 7 we prove the existence of 
these orbits. Moreover, for equal masses, we provide their complete 
classification in Section 8, for k > 0, and Section 9, for k < 0. 

Finally, in Section 10, we discuss the existence of hyperbolic homo- 
graphic solutions, which occur only for negative curvature. We prove 
that when the bodies are on the same hyperbolically rotating geodesic, 
a class of solutions we call hyperbolic Eulerian, every orbit is a hy- 
perbolic Eulerian relative equilibrium. Therefore hyperbolic Eulerian 
relative equilibria are unstable, a fact that makes them unlikely ob- 
servable candidates in a (hypothetically) hyperbolic physical universe. 



2. Equations of motion 

We consider the equations of motion on 2-dimensional manifolds of 
constant curvature, namely spheres embedded in M.^, for k > 0, and 
hyperboloid^ embedded in the Minkovski space M^, for k < 0. 



The hyperboloid corresponds to Weierstrass's model of hyperbolic geometry 
(see Appendix in [2]). 
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Consider the masses mi, m2, rris > in M'^, for k > 0, and in M^, for 
K < 0, whose positions are given by the vectors = {xi,yi,Zi), i = 
1,2,3. Let q = (qi,q2,q3) be the configuration of the system, and 
P = (Pi;P2,P3), with Pi = rriiCii, representing the momentum. We 
define the gradient operator with respect to the vector q^ as 

where a is the signature function, 



(1) a 



+1, for K> 
-1, for K < 0, 



and let V denote the operator (Vq^, Vqj, Vqg). For the 3-dimensional 
vectors a = {ax, ay, a^) and b = {bx, by, b^), we define the inner product 

(2) a b := {aj)^ + ayby + aaj)^) 
and the cross product 

(3) a b := {ayb^ - aj)y, aj)^ - aj)-^, a{a.j,by - ayb^)). 

The Hamihonian function of the system describing the motion of the 
3-body problem in spaces of constant curvature is 

H^{q, p) = r^(q, p) - U^{q), 

where 

1 ^ 

Tk{(1, p) = ^^rn^\Pi Pi)(/tqi © Qi) 

i=l 

defines the kinetic energy and 

(A) TJ (n) = V m,mj|/t|^/^/tq,0q,- 

^ ' ^^^^ ,<^<3 q.)('^q, q,) - ^(«^q. q,)^]^/' 

is the force function, —U^ representing the potential energ j13 Then the 
Hamiltonian form of the equations of motion is given by the system 



(5) 



Pi = Vq^t/«(q) - m. ^K{pi Pi)qi, i = 1,2,3, n ^ 0, 



In [2] , we showed how this expression of Uk follows from the cotangent potential 
for K 0, and that Uq is the Newtonian potential of the Euclidean problem, 
obtained as k — 0. 
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where the gradient of the force function has the expression 

(6) 
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Vq. [a(«:q, q,)(«:q, q,) - a(«q, q,-)^]^/^ 

The motion is confined to the surface of nonzero constant curvature k, 
i.e. (q, p) G T*(M^)^, where T*(M^)^ is the cotangent bundle of the 
configuration space (M^)^, and 

= {{x, y, z) eW^ \ k{x^ + + (Tz^) = 1}. 

In particular, = is the 2-dimensional sphere, and M^^^ = is 
the 2-dimensional hyperbolic plane, represented by the upper sheet of 
the hyperboloid of two sheets (see the Appendix of [2] for more details). 
We will also denote by for k > 0, and by for k <0. 

Notice that the 3 constraints given by Kq^ qj = 1, z = 1, 2, 3, imply 
that qi Pi = 0, so the 18-dimensional system ([s]) has 6 constraints. 
The Hamiltonian function provides the integral of energy, 

^«(q,p) = h, 

where h is the energy constant. Equations ([s]) also have the integrals 
of the angular momentum, 

3 

(7) ^ qi Pi = c, 

i=l 

where c = (a, /3, 7) is a constant vector. Unlike in the Euclidean case, 
there are no integrals of the center of mass and linear momentum. 
Their absence complicates the study of the problem since many of the 
standard methods don't apply anymore. 

Using the fact that /tqj q^ = 1 for i = 1, 2, 3, we can write system 
(§ as 



qj qj)[iKqi qi)qj - (nqi qj)q/ 



[a - a{Kqi 

which is the form of the equations of motion we will use in this paper. 

3. Local existence and uniqueness of Lagrangian 

solutions 



In this section we define the Lagrangian solutions of the curved 3- 
body problem, which form a particular class of homographic orbits. 
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Then, for equal masses and suitable initial conditions, we prove their 
local existence and uniqueness. 

Definition 1. A solution of equations ([S]) is called Lagrangian if, at 
every time t, the masses form an equilateral triangle that is orthogonal 
to the z axis. 

According to Definition [T| the size of a Lagrangian solution can vary, 
but its shape is always the same. Moreover, all masses have the same 
coordinate z{t), which may also vary in time, though the triangle is 
always perpendicular to the z axis. 

We can represent a Lagrangian solution of the curved 3-body problem 
in the form 

(9) q = (qi,q2,q3), with = (xi,?/^,^^), z = 1,2,3, 

xi=rcosa;, yi = rsmu, zi = z, 

a;2 = r cos(a; + 27r/3), ?/2 = ''^ sin(c<; + 27r/3), Z2 = z, 
= r cos(a; + 47r/3), = r sin(a; + 47r/3), Z3 = z, 

where z = z{t) satisfies z"^ = an^^ — crr^; a is the signature function 
defined in ([T]); r := r(t) is the size function; and tu := uj{t) is the 
angular function. 

Indeed, for every time t, we have that xf{t)+yf{t)+azf(t) = i = 
1, 2, 3, which means that the bodies stay on the surface M^, each body 
has the same z coordinate, i.e. the plane of the triangle is orthogonal 
to the z axis, and the angles between any two bodies, seen from the 
geometric center of the triangle, are always the same, so the triangle 
remains equilateral. Therefore representation (|9| of the Lagrangian 
orbits agrees with Definition [T] 

Definition 2. A Lagrangian solution of equations ([s]) is called La- 
grangian homothetic if the equilateral triangle expands or contracts, 
but does not rotate around the z axis. 

In terms of representation ([o]), a Lagrangian solution is Lagrangian 
homothetic if uj{t) is constant, but r(t) is not constant. Such orbits 
occur, for instance, when three bodies of equal masses lying initially 
in the same open hemisphere are released with zero velocities from an 
equilateral configuration, to end up in a triple collision. 

Definition 3. A Lagrangian solution of equations ([S]) is called a La- 
grangian relative equilibrium if the triangle rotates around the z axis 
without expanding or contracting. 
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In terms of representation ([o]), a Lagrangian relative equilibrium 
occurs when r(t) is constant, but u(t) is not constant. Of course, 
Lagrangian homothetic solutions and Lagrangian relative equilibria, 
whose existence we proved in p] , are particular Lagrangian orbits, but 
we expect that the Lagrangian orbits are not reduced to them. We now 
show this by proving the local existence and uniqueness of Lagrangian 
solutions that are neither Lagrangian homothetic, nor Lagrangian rel- 
ative equilibria. 

Theorem 1. In the curved 3-body problem of equal masses, for every 
set of initial conditions belonging to a certain class, the local existence 
and uniqueness of a Lagrangian solution, which is neither Lagrangian 
homothetic nor a Lagrangian relative equilibrium, is assured. 



Proof. We will check to see if equations ([8j) admit solutions of the form 
^ that start in the region z > and for which both r(t) and cu{t) are 
not constant. We compute then that 

/^Qi qj = 1 — 3Kr^/2 for i,j = 1,2, 3, with i =^ j, 



Xi = r cos u — ruj sin w, Vi = r sin u + ru cos u 

X2 = r cos (cu + — ] — rco sin ( cu + 



3 / V 3 



y2 = r sin (^u + ^ j + ru cos (^u + , 



/ 47r\ / 47r\ 

Xs = r cos yu + — j — ru sm + — j , 

/ 47r\ / 47r\ 

ys = r sm + — j +rucos [u + — j, 

(10) zi = Z2 = Zs = —arr{aKr^ — ar'^y^^'^, 

• 2 

9 9 /xT^ 

KC{i Q c{i = Kr 6j +- for i = 1,2,3, 

1 — Kr-^ 

Xi = (f — ruP') cosuj — {ruj + 2ruj) sino;, 
iji = (f — rw^) sinw + (rid + 2r6j) cosw, 

X2 = [r — ru ) cos yu + — j — [ru + 2ru) sm + — 

, 2\ 27r\ , , / 27r 

y2 = [r — ru ) sm yu + — j + [ru + 2ru) cos + — 

, yru + 2fu) sin (u+ . 

:7r\ / 47i 

— j + (rw + 2ru) cos (^w + — 



X3 = (r — rw j cos I u + 

yj, = [r — ru ) sin (^o; + — j + (rw + 2f(i') cos \ u + 
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Zi = Z2 = = —arr[aK ^ — crr^) ^^'^ — k ^r'^[aK ^ — ar"^) 



Substituting these expressions into system ([8j), we are led to the sys- 
tem below, where the double-dot terms on the left indicate to which 
differential equation each algebraic equation corresponds: 

Xi : Acosu — B smoj = 0, 

X2 : A cos (^u + ^ j - B sin (^u + ^ j = 0, 

X3 : A cos (^LO + ^ j - B sin (u + ^ j = 0, 
jji : Asinw + -B cosw = 0, 

i/2 : ^sin (^u + ^ j + B cos (^u + ^ j = 0, 

m ■ A sin (u + ^ j + B cos (u + ^ j = 0, 

^1,^2, ^3: A = 0, 

where 



A := A{t) = f - r(l - Kr^)^^ + 



Krr 24m(l — nr ) 



K,r 



2 r2(12-9Kr2)3/2' 



5 := B{t) =riu + 2ruj. 

Obviously, the above system has solutions if and only if A = B = 0, 
which means that the local existence and uniqueness of Lagrangian 
orbits with equal masses is equivalent to the existence of solutions of 
the system of differential equations 



V 

2uw 



(11) {w = - 

V = r(l — Kr"')w 



2\^,,2 _ Kru'^ 24m{l-Kr^) 

1-Kr2 r2(12-9Kr2)3/2 5 



with initial conditions r(0) = ro,w{0) = wo,i^{0) = uq, where w = 
uj. The functions r, w, and w are analytic, and as long as the initial 
conditions satisfy the conditions ro > for all k, as well as ro < 
K~^/2 for K > 0, standard results of the theory of differential equations 
guarantee the local existence and uniqueness of a solution (r, w, p) of 



equations (11), and therefore the local existence and uniqueness of a 
Lagrangian orbit with r{t) and uj{t) not constant. The proof is now 
complete. □ 
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4. Classification of Lagrangian solutions for k> 

We can now state and prove the following result: 

Theorem 2. In the curved 3-body problem with equal masses and k > 
there are five classes of Lagrangian solutions: 

(i) Lagrangian homothetic orbits that begin or end in total collision 
in finite time; 

(a) Lagrangian relative equilibria that move on a circle; 
(Hi) Lagrangian periodic orbits that are neither Lagrangian homo- 
thetic nor Lagrangian relative equilibria; 

(iv) Lagrangian non-periodic, non-collision orbits that eject at time 
—oo, with zero velocity, from the equator, reach a maximum distance 
from the equator, which depends on the initial conditions, and return 
to the equator, with zero velocity, at time +oo. 

None of the above orbits can cross the equator, defined as the great 
circle of the sphere orthogonal to the z axis. 

(v) Lagrangian equilibrium points, when the three equal masses are 
fixed on the equator at the vertices of an equilateral triangle. 

The rest of this section is dedicated to the proof of this theorem. 



Let us start by noticing that the first two equations of system (11) 
imply that w = — which leads to 

c 

where c is a constant. The case c = can occur only when w = 0, which 
means u = 0. Under these circumstances the angular velocity is zero, so 
the motion is homothetic. These are the orbits whose existence is stated 
in Theorem[2](i). They occur only when the angular momentum is zero, 
and lead to a triple collision in the future or in the past, depending on 
the sense of the velocity vectors. 



For the rest of this section, we assume that c 7^ 0. Then system (11 ) 
takes the form 



(12) 



_ c^(l-Kr^) _ Kru'^ 24m(l-Kr^) 



Notice that the term of the last equation arises from the deriva- 



tives zi,Z2,Z'i in (10). But these derivatives would be zero if the equi- 
lateral triangle rotates along the equator, because r is constant in this 

2 

case, so the term -fz^ vanishes. Therefore the existence of equilateral 
relative equilibria on the equator (included in statement (ii) above), 
and the existence of equilibrium points (stated in (v)) — results proved 



10 



Florin Diacu and Ernesto Perez-Chavela 



in [2] — are in agreement with the above equations. Nevertheless, the 

2 

term stops any orbit from crossing the equator, a fact mentioned 



1— /tr^ 

before statement fv) of Theorem |2 



Understanding system (12) is the key to proving Theorem |2| We 
start with the following facts: 

Lemma 1. Assume n^m > Q and c 7^ 0. Then for h^I'^c? — {9>/\/?))m < 
0, system (12) has two fixed points, while for k^/^c^ — {8/^/3)m > it 
has one fixed point. 



Proof. The fixed points of system (12) are given by f = = z>. Substi- 
tuting u = in the second equation of (12), we obtain 

0. 



1 — Kr^ fc^ 



24m 



(12 -9Kr2)3/2 



The above remarks show that, for k > 0, r = ' is a fixed point, 
which physically represents an equilateral relative equilibrium moving 
along the equator. Other potential fixed points of system ( 12 ) are given 
by the equation 

c^(12 - dKr^/^ = 24mr, 
whose solutions are the roots of the polynomial 

(13) 729c^/tV^ - 2916c^kV + 144(27c^k + im^y - 1728. 
Writing x = r"^ and assuming k, > 0, this polynomial takes the form 

(14) p{x) = 729k^x^ - 2916cV^x2 + 144(27cV + 4m^)a; - 1728, 
and its derivative is given by 

(15) p'{x) = 2187c^K^x^ - 5832cV2x + 144(27c^k + 4m^). 

The discriminant of p' is — 5038848c^fi;^m^ < 0. 

By Descartes's rule of signs, p can have one or three positive roots. 
If p has three positive roots, then p' must have two positive roots, but 
this is not possible because its discriminant is negative. Consequently 
p has exactly one positive root. 

For the point (r, = (rg, 0) to be a fixed point of equations (12), tq 
must satify the inequalities < tq < If we denote 

(16) gir) 



c 
r 



24m 



(12-9ft;r2)3/2' 
we see that, for k > 0, (7 is a decreasing function since 



(17) 



d_ 
dr 



9[r) 



648mKr 



;i2 - 9«:r2)5/2 



< 0. 
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When r — 0, we obviously have that g{r) > since we assumed 
c 7^ 0. When r — > we have g{r) — )■ k-'^/^c^ — (8/v^)m. If 

— (8/-\/3)'Ti > 0, then tq > so (ro,0) is not a fixed point. 

Therefore, assuming c 7^ 0, a necessary condition that (rg, 0) is a fixed 
point of system (12) with < tq < k"^/^ is that 



K^/^c^ - {8/V3)m < 0. 

For K-'^/^c^ — (8/-\/3)'Ti > 0, the only fixed point of system (12) is (r, u) = 
(k~^/^,0). This conclusion completes the proof of the lemma. □ 

4.1. The flow in the (r, u) plane for k > 0. We will now study the 
fiow of system (12) in the (r, u) plane for k > 0. At every point with 



z/ 7^ 0, the slope of the vector field is given by i.e. by the ratio 



h(r, u), where 
h{r, v) = 



Kr 



24m (1 



KT 



z/r2(12- W)3/2- 



Since h{r,—u) = —h(r,iy), the fiow of system (12) is symmetric with 
respect to the r axis for r G (0,^"^/^]. Also notice that, except for 
the fixed point {k,~^^'^,0), system (12) is undefined on the lines r = 
and r = k~^^^. Therefore the fiow of system (12) exists only for points 
(r, z/) in the band (0, k"^/^) x M and for the point {k~^/'^, 0). 

Since r = z/, no interval on the r axis can be an invariant set for 



system (12). Then the symmetry of the fiow relative to the r axis 



implies that orbits cross the r axis perpendicularly. But since g{r) 7^ 
at every non- fixed point, the fiow crosses the r axis perpendicularly 
everywhere, except at the fixed points. 

Let us further treat the case of one fixed point and the case of two 
fixed points separately. 

4.1.1. The case of one fixed point. A single fixed point, namely 
(k~-^/2,0), appears when k'/^c^ - (8 /V3)m > 0. Then the function g, 
which is decreasing, has no zeroes for r G (0, k"^/^), therefore g{r) > 
in this interval, so the fiow always crosses the r axis upwards. 



For z^ 7^ 0, the right hand side of the second equation of (12) can be 
written as 



(18) 
where 

(19) 



Gir, v) = gi{r)g{r) + g2{r, u), 



Kr 



and g2{r,u) 



1 — nr"^ 
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But £gi{r) = -2/r3 < and |(72(r,i/) = -^^^^ < 0. So, like 

g, the functions gi and g2 are decreasing in (0, with gi,g > 0, 

therefore G is a decreasing function as weU. Consequently, for u = 
constant > 0, the slope of the vector field decreases from +oo at r = 
to — oo ar r = k"^/^. For v = constant < 0, the slope of the vector 
field increases from — oo at r = to +oo at r = 

This behavior of the vector field forces every orbit to eject downwards 
from the fixed point, at time t = —oo and with zero velocity, on a 
trajectory tangent to the line r = k^^/^, reach slope zero at some 
moment in time, then cross the r axis perpendicularly upwards and 
symmetrically return with final zero velocity, at time t = +oo, to 
the fixed point (see Figure [l]^a)). So the fiow of system (12) consists 
in this case solely of homoclinic orbits to the fixed point (k~^/^,0), 
orbits whose existence is claimed in Theorem [2] (iv). Some of these 
trajectories may come very close to a total collapse, which they will 
never reach because only solutions with zero angular momentum (like 
the homothetic orbits) encounter total collisions, as proved in [3]. 

So the orbits cannot reach any singularity of the line r = 0, and 
neither can they begin or end in a singularity of the line r = 
The reason for the latter is that such points are or the form (k^^/^, u) 
with u ^ 0, therefore r 7^ at such points. But the vector field 
tends to infinity when approaching the line r = so the fiow 

must be tangent to it, consequently r must tend to zero, which is a 
contradiction. Therefore only homoclinic orbits exist in this case. 




Figure 1. A sketch of the fiow of system ([12]) for (a) 
K = c = l,m = 0.24, typical for one fixed point, and (b) 
K = c=l,m = 4, typical for two fixed points. 
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4.1.2. The case of two fixed points. Two fixed points, 
and (ro,0), with < tq < k^^^^, occur when k^/^c^ — {8/\/3)m < 0. 
Since g is decreasing in the interval (0, we can conclude that 

g{t) > for t e (0,ro) and g{t) < for t G (ro,/t"^/^). Therefore 



the flow of system (12) crosses the r axis upwards when r < Tq, but 
downwards for r > rg (see Figure [l](b)). 

The function G{r,u), defined in ( jlsj ), fails to be decreasing in the 
interval (0, along lines of constant u, but it has no singularities 

in this interval and still maintains the properties 

lim G(r, u) = +oo and lim G(r, v) = — oo. 

Therefore G must vanish at some point, so due to the symmetry of 
the vector field with respect to the r axis, the fixed point {tq, 0) is 
surrounded by periodic orbits. The points where G vanishes are given 
by the nullcline z> = 0, which has the expression 

2 (1 - nr'^y 24m 

Kr^ r (12 — 9Kr^)^/^ 

This nullcline is a disconnected set, formed by the fixed point (/t^^/^, 0) 
and a continuous curve, symmetric with respect to the r axis. Indeed, 
since the equation of the nullcline can be written as u'^ = ''"^"^'^3 g{r), 
and lim^^(^-i/2)- g{r) = /t^/^c^ — (8/-\/3)m < in the case of two fixed 
points (as shown in the proof of Lemma [l|, only the point (k^^/^,0) 
satisfies the nullcline equation away from the fixed point (ro, 0). 

The asymptotic behavior of G near r = also forces the flow 

to produce homoclinic orbits for the fixed point {k~^^'^,0), as in the 
case discussed in Subsection 4.1.1. The existence of these two kinds 
of solutions is stated in Theorem [2] (iii) and (iv), respectively. The 
fact that orbits cannot begin or end at any of the singularities of the 
lines r = or r = k"^/^ follows as in Subsection 4.1.1. This remark 
completes the proof of Theorem 2. 



5. Classification of Lagrangian solutions for k < 

We can now state and prove the following result: 

Theorem 3. In the curved 3-body problem with equal masses and k, < 
there are eight classes of Lagrangian solutions: 

(i) Lagrangian homothetic orbits that begin or end in total collision 
in finite time; 

(a) Lagrangian relative equilibria, for which the bodies move on a 
circle parallel with the xy plane; 
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(Hi) Lagrangian periodic orbits that are not Lagrangian relative equi- 
libria; 

(iv) Lagrangian orbits that eject at time — oo from a certain rela- 
tive equilibrium solution s (whose existence and position depend on the 
values of the parameters) and returns to it at time +00; 

(v) Lagrangian orbits that come from infinity at time —00 and reach 
the relative equilibrium s at time +00; 

(vi) Lagrangian orbits that eject from the relative equilibrium s at 
time —00 and reach infinity at time +00; 

(vii) Lagrangian orbits that come from infinity at time —00 and sym- 
metrically return to infinity at time +00, never able to reach the La- 
grangian relative equilibrium s; 

(via) Lagrangian orbits that come from infinity at time —00, reach a 
position close to a total collision, and symmetrically return to infinity 
at time +00. 

The rest of this section is dedicated to the proof of this theorem. No- 
tice first that the orbits described in Theorem |3] (i) occur for zero angu- 
lar momentum, when c = 0, as for instance when the three equal masses 
are released with zero velocities from the Lagrangian configuration, a 
case in which a total collapse takes place at the point (0,0, 
Depending on the initial conditions, the motion can be bounded or un- 
bounded. The existence of the orbits described in Theorem |3] (ii) was 
proved in [2]. To address the other points of Theorem [3| and show that 
no other orbits than the ones stated there exist, we need to study the 



flow of system (12) for k < 0. Let us first prove the following fact: 



Lemma 2. Assume k < 0,m > 0, and c 7^ 0. Then system (12) has 
no fixed points when 27c^k, + 4m^ < 0, and can have two, one, or no 
fixed points when 27c'*k + 4m^ > 0. 



Proof. The number of fixed points of system (12) is the same as the 
number of positive zeroes of the polynomial p defined in ( 14 ) . If 27c^k-|- 
4m^ < 0, all coefficients of p are negative, so by Descartes's rule of signs, 
p has no positive roots. 

Now assume that 27c'*/«+4m^ > 0. Then the zeroes of p are the same 
as the zeroes of the monic polynomial (i.e. with leading coefficient 1): 



p[x) =x^ - Ak'^x^ + [48k~^ + (64/81)c~^K-W]x - (64/27)/t-^ 



obtained when dividing p by the leading coefficient. But a monic cubic 
polynomial can be written as 

— (oi + 02 + 03)0;^ + (aia2 + 0203 + a2,ai)x — 010203, 
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where ai, 02, and are its roots. One of these roots is always real and 
has the opposite sign of —010203. Since the free term of p is positive, 
one of its roots is always negative, independently of the allowed values 
of the coefficients k, m, c. Consequently p can have two positive roots 
(including the possibility of a double positive root) or no positive root 



at all. Therefore system (12) can have two, one, or no fixed points. As 



we will see later, all three cases occur. □ 



We further state and prove a property, which we will use to establish 
Lemma HI 

Lemma 3. Assume k < 0, m > 0, c 7^ 0, let (r^,,0) be a fixed point 



of system (|12|), and consider the function g defined in (|16|). Then 



^g{r*) = if and only if = (-|:)^/^- Moreover, ^igir^) > 0. 



Proof. Since (r=i,,0) is a fixed point of system (12), it follows that 
(7(r*) = 0. Then it follows from relation (16) that (12 — Q/tr^"!^/^ = 



24mr^,/c^. Substituting this value of (12 — into the equation 

^g{r^) = 0, which is equivalent to 

648m/tr^, 
(12 - 9Kr2)5/2 = 

it follows that 27k/(12 - Q/tr^) = -1/rl Therefore r, = (-|:)^/^. 
Obviously, for this value of r^,, g{r^:) = 0, so the first part of Lemma |3] 
is proved. To prove the second part, substitute = (— ^)^^^ into the 
equation g{r^) = 0, which is then equivalent with the relation 

(20) 9V3c^{-kY/^ - 4m = 0. 
Notice that 

rf2 ^ 2c2 648mK 29160mK^r^ 
J9{r 



dr^""- ' r3 (12-W)5/2 (12-9fi:r2)7/2- 



Substituting for r^, = (— ^)^/^ in the above equation, and using (20) 



we are led to the conclusion that j^gijr^) = — (2\/3 + 6y/2)mK/9y/6, 
which is positive for n < 0. This completes the proof. □ 

The following result is important for understanding a qualitative 



aspect of the fiow of system (12), which we will discuss later in this 
section. 

Lemma 4. Assume k < 0, m > 0, c 7^ 0, and let (r^,, 0) be a fixed point 
|-G(r„0) = 0, then 



of system (Il2]). // |-G'(r^,0) = 0, then ^G(r^,0) > 0, where G is 



defined in (18). 



16 



Florin Diacu and Ernesto Perez-Chavela 



Proof. Since (r*, 0) is a fixed point of (12), G{r^:,0) = 0. But for /t < 



we have gi{r^) > 0, so necessarily g{r^:) = 0. Moreover, -^giir^) ^ 0, 
and since ^5-2 (r, u) = - ^^j^z^^Syr- , it follows that §:g2{r*, 0) = 0. But 
do d d d 

so the condition ^G'(r^,,0) = implies that -^g{r^) = 0. By Lemma 
, n = (— J:)"'^''^ and ^g{r^) > 0. Using now the fact that 

it follows that ^G{r^, 0) = g\{r^)-^g{r^). Since Lemma [s] implies that 
^g{r^ > 0, and we know that > 0, it follows that ^G(r^,, 0) > 

0, a conclusion that completes the proof. □ 

5.1. The flow in the (r, u) plane for k, < 0. We will now study the 



flow of system (12) in the (r, plane for k < 0. As in the case /t > 0, 
and for the same reason, the flow is symmetric with respect to the 
r axis, which it crosses perpendicularly at every non-fixed point with 
r > 0. Since we can have two, one, or no fixed points, we will treat 
each case separately. 

5.1.1. The case of no fixed points. No fixed points occur when g{r) 
has no zeroes. Since g{r) — >■ 00 as r — )■ with r > 0, it follows that 
g{r) > 0. Since gi{r) and g2{r,p) are also positive, it follows that 
G{r,v) > for r > 0. But h{r,v) = G{r,v)/v. Then h{r,v) > for 
z/ > and h{r, u) < for u < 0, so the flow comes from infinity at time 
—00, crosses the r axis perpendicularly upwards, and symmetrically 
reaches infinity at time +00 (see Figure ga)). These are orbits as in 
the statements of Theorem |3] (vii) and (viii) but without any reference 
to the Lagrangian relative equilibrium s. 

5.1.2. The case of two fixed points. In this case, the function g 



defined at ( 16 ) has two distinct zeroes, one for r = ri and the other for 
r = r2, with < ri < r2. In Theorem [3} we denoted the fixed point 
(r2,0) by s. Moreover, g{r) > for r G (0,ri) U (r2,oo), and g{r) < 
for r G (ri, r2). Therefore the vector field crosses the r axis downwards 
between ri and r2, but upwards for r < ri as well as for r > r2. 

To determine the behavior of the flow near the fixed point (ri,0). 



we linearize system (12). For this let F(r,iy) = u to he the right 



hand side of the first equation in (12), and notice that ^(ri,0) = 1, 



|^(ri,0) = 0, and ^(ri,0) = 0. Since, along the r axis, G{r,0) 
is positive for r < ri, but negative for ri < r < r2, it follows that 
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kft, }f //// y,/^^ 

y. , ^ ' . h ' . . . . ,. f ...^ ; . 



^1*1 ^ftooH "401)0 '\ ^Wdt '^"^ ^iio\v V lodoo 



ftp -v^^ 



\vi \m\'' \\\\\ \^ -^ix^^ 




Figure 2. A sketch of the flow of system (12) for (a) 
K = — 2,c = 1/3, and m = 1/2, typical for no fixed 
points; (b) k = —0.3, c = 0.23, and m = 0.12, typical for 
two fixed points, which are in this case on the line u = 
at approximately ri = 1.0882233 and = 2.0007055. 



either ^(ri,0) < or ^(ri,0) = 0. But according to Lemma |4| if 

^(ri,0) = 0, then |^G(r^,0) > 0, so G{r,0) is convex up at (ri,0). 
Then G{r, 0) cannot not change sign when r passes through ri along 
the line z/ = 0, so the only existing possibility is 0) < 0. 

The eigenvalues of the linearized system corresponding to the fixed 
point (ri, 0) are then given by the equation 



(21) 



Since ^ 

or 



det 



dr 



-A 1 
ri,0) -A 



0. 



(ri,0) is negative, the eigenvalues are purely imaginary, so 



(ri,0) is not a hyperbolic fixed point for equations (12). Therefore 
this fixed point could be a spiral sink, a spiral source, or a center for 



the nonlinear system. But the symmetry of the flow of system (12) 



with respect to the r axis, and the fact that, near ri, the flow crosses 
the r axis upwards to the left of ri, and downwards to the right of ri, 
eliminates the possibility of spiral behavior, so (ri,0) is a center (see 
Figure gb)). 

We can understand the generic behavior of the flow near the isolated 
fixed point (r2, 0) through linearization as well. For this purpose, notice 
that |^(r2,0) = 1, |^(r2,0) = 0, and g(r2,0) = 0. Since, along the r 
axis, G{r) is negative for ri < r < r2, but positive for r > r2, it follows 
that ^(r2,0) > or ^(r2,0) = 0. But using Lemma W] the same way 
we did above for the fixed point (ri, 0), we can conclude that the only 
possibility is ^(r2,0) > 0. 
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The eigenvalues corresponding to the fixed point (r2, 0) are given by 
e eauation 



the equation 

(22) det 



dG, 



0. 



Consequently the fixed point (r2,0) is hyperbohc, its two eigenvalues 
are Ai > and A2 < 0, so (r2, 0) is a saddle. 

Indeed, for small z/ > 0, the slope of the vector field decreases to —00 
on lines r = constant, with ri < r < r2, when u tends to 0. On the 
same lines, with r > r2, the slope decreases from +00 as u increases. 
This behavior gives us an approximate idea of how the eigenvectors 
corresponding to the eigenvalues Ai and A2 are positioned in the ru 
plane. 

On lines of the form u = rjr, with rj > 0, the slope h{r, v) of the 
vector field becomes 



h{r, rjr) 



K.2 



^^3 



KV (? 24m (1 — nr^ 



:i2-9ft:r2)3/2j 1 



KT 



2 



So, as r tends to 00, the slope h{r,rir) tends to rj. Consequently the 
vector field doesn't bound the fiow with negative slopes, and thus allows 
it to go to infinity. 

With the fixed point (ri,0) as a center, the fixed point (^2,0) as a 
saddle, and a vector field that doesn't bound the orbits as r — )■ 00, the 
fiow must behave qualitatively as in Figure [2|^b). 

This behavior of the fiow proves the existence of the following types 
of solutions: 

(a) periodic orbits around the fixed point (ri,0), corresponding to 
Theorem |3] (iii); 

(b) a homoclinic orbit to the fixed point (ri,0), corresponding to 
Theorem |3] (iv); 

(c) an orbit that tends to the fixed point (r2,0), corresponding to 
Theorem |3] (v) ; 

(d) an orbit that is ejected from the fixed point (r2, 0), corresponding 
to Theorem [3] (vi); 

(e) orbits that come from infinity in the direction of the stable man- 
ifold of (r2, 0) and inside it, hit the r axis to the right of r2, and return 
symmetrically to infinity in the direction of the unstable manifold of 
(r2, 0); these orbits correspond to Theorem |3] (vii); 

(f) orbits that come from infinity in the direction of the stable man- 
ifold of (r'2,0) and outside it, turn around the homoclinic loop, and 
return symmetrically to infinity in the direction of the unstable mani- 
fold of (r2, 0); these orbits correspond to Theorem [s] (viii). 
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Since no other orbits show up, the proof of this case is complete. 



5.1.3. The case of one fixed point. We left the case of one fixed 
point at the end because it is non-generic. It occurs when the two 
fixed points of the previous case overlap. Let us denote this fixed point 
by (t'0,0). Then the function g{r) is positive everywhere except at the 
fixed point, where it is zero. So near rg, g is decreasing for r < tq 
and increasing for r > tq, and the r axis is tangent to the graph of 
g. Consequently, ^(ro,0) = 0, and the eigenvalues obtained from 



equation (21) are Ai = A2 = 0. In this degenerate case, the orbits near 
the fixed point influence the asymptotic behavior of the flow at (ro, 0). 
Since the flow away from the fixed point looks very much like in the 
case of no fixed points, the only difference between the flow sketched 
in Figure |2]^a) and the current case is that at least an orbit ends at 
(ro,0), and at least another orbit one ejects from it. These orbits are 
described in Theorem [sl fiv) and (v). 



The proof of Theorem |3] is now complete. 



6. Mass equality of Lagrangian solutions 

In this section we show that all Lagrangian solutions that satisfy 
Definition [1] must have equal masses. In other words, we will prove the 
following result: 

Theorem 4. In the curved 3-body problem, the bodies of masses mi, m2, 
7713 can lead to a Lagrangian solution if and only if rrii = m2 = m^. 

Proof. The fact that three bodies of equal masses can lead to La- 
grangian solutions for suitable initial conditions was proved in The- 
orem [ij So we will further prove that Lagrangian solutions can occur 
only if the masses are equal. Since the case of relative equilibria was 
settled in [2], we need to consider only the Lagrangian orbits that are 
not relative equilibria. This means we can treat only the case when 
r(t) is not constant. 

Assume now that the masses are mi, 7712,^3, and substitute a solu- 
tion of the form 



Xi = r cosu, 

X2 = r cos(a; + 27r/3), 

X3 = r cos{ijJ + 47r/3), 



yi = r sincu, Zi 
?/2 = r sin(cj + 27r/3), Z2 
?/3 = r sin(a; + 47r/3), Z3 
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into the equations of motion. Computations and a reasoning similar to 
the ones performed in the proof of Theorem [T] lead us to the system: 

9, 9 Krr^ 12(mi + molfl — Kr^) 
r-r{l- kt')0? + ^ + \' l\~,zi2 = 0' 

9, 9 Krr"^ 12(m2 + mq)(l — /tr^) 
r-r{l- Kr')u' + + \' ' = 0, 

9 9 K,rP 12(m3 + mi)(l — /tr^) 
f — r(l — Kr )Cj + r + " 



l-/s:r2 r2(12 - 9Kr2)3/2 

.. 4y3(mi -ms) 

+ 2ruj - —— — , = 0, 

.. , o - • 4y3(m2 -ms) 
+ 2ruj - —— — ^.3,^ = 0, 

.. 4y3(m3 - mi) 

+ 2rw - —— — ,^ = 0, 

which, obviously, can have solutions only if mi = m2 = m^. This 
conclusion completes the proof. □ 



7. Local existence and uniqueness of Eulerian solutions 

In this section we define the Eulerian solutions of the curved 3-body 
problem and prove their local existence for suitable initial conditions 
in the case of equal masses. 

Definition 4. A solution of equations ([s]) is called Eulerian if, at ev- 
ery time instant, the bodies are on a geodesic that contains the point 
(0,0,|«:|-V2|). 

According to Definition |4| the size of an Eulerian solution may 
change, but the particles are always on a (possibly rotating) geodesic. 
If the masses are equal, it is natural to assume that one body lies at 
the point (0,0, |k|~^/2|)^ while the other two bodies find themselves at 
diametrically opposed points of a circle. Thus, in the case of equal 
masses, which we further consider, we ask that the moving bodies have 
the same coordinate 2, which may vary in time. 

We can thus represent such an Eulerian solution of the curved 3-body 
problem in the form 



(23) q = (qi,q2,q3), with = (0:^,1/^,2;^), i = 1,2,3, 
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xi = 0, yi = 0, zi = (ctk) 

X2 = rcosu, y2 = rsmco, Z2 = z, 

X3 = —rcosu, ?/3 = — rsina;, = z, 

where z = z{t) satisfies = an^^ — ar"^ = {aK,)~^{l — nr"^); a is the 
signature function defined in ([T|; r := r{t) is the size function; and 
u := u(t) is the angular function. 

Notice that, for every time t, we have x^(t)+yf(t)+azf{t) = i = 
1, 2, 3, which means that the bodies stay on the surface M^. Equations 



(23) also express the fact that the bodies are on the same (possibly 



rotating) geodesic. Therefore representation (23) of the Eulerian orbits 



agrees with Definition |4] in the case of equal masses. 



Definition 5. An Eulerian solution of equations ([8j) is called Euler- 
ian homothetic if the configuration expands or contracts, but does not 
rotate. 



In terms of representation (23), an Eulerian homothetic orbit for 
equal masses occurs when w(t) is constant, but r(t) is not constant. If, 
for instance, all three bodies are initially in the same open hemisphere, 
while the two moving bodies have the same mass and the same z co- 
ordinate, and are released with zero initial velocities, then we are led 
to an Eulerian homothetic orbit that ends in a triple collision. 



Definition 6. An Eulerian solution of equations (|8j) is called an Euler- 
ian relative equilibrium if the configuration of the system rotates without 
expanding or contracting. 



In terms of representation (23 ), an Eulerian relative equilibrium orbit 
occurs when r(t) is constant, but oj(t) is not constant. Of course, 
Eulerian homothetic solutions and elliptic Eulerian relative equilibria, 
whose existence we proved in [2] , are particular Eulerian orbits, but we 
expect that the Eulerian orbits are not reduced to them. We now show 
this fact by proving the local existence and uniqueness of Eulerian 
solutions that are neither Eulerian homothetic, nor Eulerian relative 
equilibria. 

Theorem 5. In the curved 3-body problem of equal masses, for every 
set of initial conditions belonging to a certain class, the local existence 
and uniqueness of an Eulerian solution, which is neither homothetic 
nor a relative equilibrium, is assured. 



Proof. To check whether equations (\8^ admit solutions of the form (23) 



that start in the region z > and for which both r{t) and uj{t) are not 
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constant, we first compute that 

Kqi q2 = KQi qa = (1 - /tr^)^/^ 
qs = 1 - 
Xi = 0, yi = 0, 

X2 = r cos cu — rcu sin cu, ?/2 = r sin a; + rcj cos cu, 

= — f cos + rw sin u, y2 = —r sin uj — ruj cos w, 

crrr 



ii = 0, 



Kqa qa 



2^2 = 2:3 = 
Kqi qi = 0, 
Kqs qa 



2 • 2 I 



2 ' 



a;2 = [r 
i/2 = (r — ruj 
X3 = —[r — ru 



xi = jji = zi = 0, 
2^ 



rw^j cosw — {ru + 27^^) sin a; 
sin a; + (rio + 2r6j) cosu. 

2 



— if — ru'^] 



cosu + {ruj + 2fuj) sinw, 
sinw — {ruj + 2ruj) cosu, 



Z2 



or~] ^'~-K ^f'^ioK ^-ar^) 



Substituting these expressions into equations ([8j), we are led to the 
system below, where the double-dot terms on the left indicate to which 
differential equation each algebraic equation corresponds: 

C cos UJ — D sin a; = 0, 
C sinu + D cos a; = 0, 
C = 0, 



X2,a;3 

^2, ^3 



where 



(:7:=C'(t) =f-r(l 



nrf"^ 



+ 



m(5 — 4Kr^) 



1 — /tr^ ' 4r^(l — /tr^)^/^ ' 

(The equations corresponding to and zi are identities, so they 

don't show up). The above system has solutions if and only if C = 
= 0, which means that the existence of Eulerian homographic orbits 
of the curved 3-body problem with equal masses is equivalent to the 
existence of solutions of the system of differential equations: 
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with initial conditions r(0) = ro,w(0) = wq, z^(O) = uq, where w = 
uj. The functions r, w, and w are analytic, and as long as the initial 
conditions satisfy the conditions ro > for all k, as well as tq < 
for K > 0, standard results of the theory of differential equations 
guarantee the local existence and uniqueness of a solution (r, w, u) of 



equations (24), and therefore the local existence and uniqueness of 
an Eulerian orbit with r(t) and uj{t) not constant. This conclusion 
completes the proof. □ 

8. Classification of Eulerian solutions for k > 

We can now state and prove the following result: 

Theorem 6. In the curved 3-body problem with equal masses and n > 
there are three classes of Eulerian solutions: 

(i) homothetic orbits that begin or end in total collision infinite time; 

(a) relative equilibria, for which one mass is fixed at one pole of the 
sphere while the other two move on a circle parallel with the xy plane; 

(Hi) periodic homographic orbits that are not relative equilibria. 

None of the above orbits can cross the equator, defined as the great 
circle orthogonal to the z axis. 

The rest of this section is dedicated to the proof of this theorem. 
Let us start by noticing that the first two equations of system (24) 
imply that w = — which leads to 

c 

where c is a constant. The case c = can occur only when w = 0, which 
means u = 0. Under these circumstances the angular velocity is zero, 
so the motion is homothetic. The existence of these orbits is stated in 
Theorem [6] (i). They occur only when the angular momentum is zero, 
and lead to a triple collision in the future or in the past, depending 
on the direction of the velocity vectors. The existence of the orbits 
described in Theorem [6] (ii) was proved in |2]. 



For the rest of this section, we assume that c 7^ 0. System (24) is 
thus reduced to 



V 



(25) _ _ 

To address the existence of the orbits described in Theorem 161(111), and 
show that no other Eulerian orbits than those of Theorem l6| exist for 



ft > 0, we need to study the flow of system (25) for k > 0. Let us first 
prove the following fact: 
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Lemma 5. Regardless of the values of the parameters m,K > 0, and 
c 7^ 0, system (25) has one fixed point (tq, 0) with < ro < 



Proof. The fixed points of system (25) are of the form (r, 0) for all 
values of r that are zeroes of u(r), where 



^ ^ ^ ^ r 4(1 - Kr2)i/2 

But finding the zeroes of u{r) is equivalent to obtaining the roots of 
the polynomial 

16K^(c^/t + m'^y^ - 8/t(6cV + 5m^)r^ + (48c^k + 25m\^ - 16c^ 
Denoting x = r^, this polynomial becomes 

q{x) = 16K^{c'^K + m^)x^ -8K{6c^K + 5m^)x^ + {48c'^K + 25m^)x-16c\ 

Since k > 0, Descarte's rule of signs implies that q can have one or 
three positive roots. The derivative of q is the polynomial 

(27) q'{x) = 48k^(c^k + m^)x^ - 1Qk{6c^k + 5m^)x + ASc^k + 2hm^ , 

whose discriminant is ^'^K^vn^i^Xc^K + 25m'^). But, as k > 0, this 
discriminant is always positive, so it offers no additional information 
on the total number of positive roots. 

To determine the exact number of positive roots, we will use the 
resultant of two polynomials. Denoting by Oj, z = 1, 2, . . . , C, the roots 
of a polynomial P, and by bj,j = 1, 2, ... those of a polynomial Q, 
the resultant of P and Q is defined by the expression 

C 5 

Res{P,Q) = l[l[{a,-bj). 

i=i j=i 

Then P and Q have a common root if and only if Res[P, Q] = 0. 
Consequently the resultant of q and q' is a polynomial in k, c, and m 
whose zeroes are the double roots of q. But 

Res(g, q') = 1024c^ K^m^(c^K + m^)(108c^K + 125m2). 

Then, for m, n > and c 7^ 0, Res[g, q'] never cancels, therefore q has 
exactly one positive root. Indeed, should q have three positive roots, a 
continuous variation of k, m, and c, would lead to some values of the 
parameters that correspond to a double root. Since double roots are 
impossible, the existence of a unique equilibrium (ro, 0) with ro > is 
proved. To conclude that ro < k^^^"^ for all m, k > and c 7^ 0, it is 
enough to notice that \imr^Qu{r) = +00 and lim^^^-1/2 ^(r) = —00. 
This conclusion completes the proof. □ 
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8.1. The flow in the (r, u) plane for k > 0. We can now study the 
flow of system (25) in the (r, z/) plane for k, > 0. The vector field is 
not defined along the lines r = and r = k^^/^, so it lies in the band 
(0, X M. Consider now the slope ^ of the vector field. This slope 

is given by the ratio ^ = f (r, i/), where 



(2 



f (r, u) 



m(5 — 4Kr^ 



1 



But V is odd with respect to z/, i.e. v{r, —v) = — t>(r, z/), so the vector 
field is symmetric with respect to the r axis. 




"mr :^ 



Figure 3. A sketch of the flow of system (25) for k 



1, c = 2, and m = 2, typical for Eulerian solutions with 
K > 0. 

Since limr_5.o '^(r) = +oo and lim^_^^-i/2 f(r) = — oo, the flow crosses 
the r axis perpendicularly upwards to the left of tq and downwards 



to its right, where (ro,0) is the flxed point of the system (25) whose 
existence and uniqueness we proved in Lemma [Sj But the right hand 
side of the second equation in ( [25) ) is of the form 

(29) W{r,iy) = u{r)/r'^ + g2{r,u), 

where §2 was deflned earlier as g2{r, u) = — j^^, while u{r) was deflned 



in (|26j). Notice that 

lim W{r, v) 



+00 and lim W{r, v) = —oo. 



Moreover, W{rQ, 0) = 0, and W has no singularities for r G (0, 
Therefore the flow that enters the region z/ > to the left of tq must 
exit it to the right of the flxed point. The symmetry with respect to the 
r axis forces all orbits to be periodic around {tq, 0) (see Figure [s]). This 
proves the existence of the solutions described in Theorem [6] (iii), and 
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shows that no orbits other than those in Theorem [6] occur for k > 0. 
The proof of Theorem |6] is now complete. 

9. Classification of Eulerian solutions for k <0 

We can now state and prove the following result: 

Theorem 7. In the curved 3-body problem with equal masses and k > 
there are four classes of Eulerian solutions: 

(i) Eulerian homothetic orbits that begin or end in total collision in 
finite time; 

(a) Eulerian relative equilibria, for which one mass is fixed at the 
vertex of the hyperboloid while the other two move on a circle parallel 
with the xy plane; 

(Hi) Eulerian periodic orbits that are not relative equilibria; the line 
connecting the two moving bodies is always parallel with the xy plane, 
but their z coordinate changes in time; 

(iv) Eulerian orbits that come from infinity at time —oo, reach a 
position when the size of the configuration is minimal, and then return 
to infinity at time +oo. 

The rest of this section is dedicated to the proof of this theorem. 

The homothetic orbits of the type stated in Theorem [7](i) occur only 
when c = 0. Then the two moving bodies collide simultaneously with 
the fixed one in the future or in the past. Depending on the initial 
conditions, the motion can be bounded or unbounded. 

The existence of the orbits stated in Theorem [T] (ii) was proved in 
[2] . To prove the existence of the solutions stated in Theorem [T] (iii) 
and (iv), and show that there are no other kinds of orbits, we start 
with the following result: 

Lemma 6. In the curved three body problem with k < 0, the polynomial 
q defined in the proof of Lemma\^ has no positive roots for c^K+m?' < 0, 
but has exactly one positive root for c^k + > 0. 

Proof. We split our analysis in three different cases depending on the 
sign of c^K + m^: 

(1) c'^K + = 0. In this case q has form %Km?x'^ + 2?)C^kx — IGc"^, 
a polynomial that does not have any positive root. 

(2) c'^K + m? < 0. Writing Qc^K + bm^ = 6(c^K + m^) — m^, we see that 
the term of q corresponding to is always negative, so by Descartes's 
rule of signs the number of positive roots depends on the sign of the 
coefficient corresponding to x, i.e. 48c^K + 25m^ = 48(c^K + m^) — 23m^, 
which is also negative, and therefore q has no positive root. 

(3) c'^K + > 0. This case leads to three subcases: 
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- if Gc'^K + 5m^ < 0, then necessarily 48c^fi; + 25m^ < and, so q has 
exactly one positive root; 

- if 6c^K, + 5m^ > and 48c^/t + 25m^ < 0, then q has one change of 
sign and therefore exactly one positive root; 

- if 48c'^K + 25m^ > 0, then all coefficients, except for the free term, 
are positive, therefore q has exactly one positive root. 

These conclusions complete the proof. □ 

The following result will be used towards understanding the case 



when system (25) has one fixed point. 



Lemma 7. Regardless of the values of the parameters k, < 0,m > 0, 



^-W{r^,0) = 0, where W is defined in p9| ) 
Proof. Since u{r^) = 0, -^g2{r^,0) = 0, and 



and c 7^ 0, there is no fixed point, (r*,0), of system (25) for which 

dr ^ 

a 



dW 
dr 



d 



d 



(r,z/) = -(2/r )u{r) + (1/r )-ru{r) + —g2{r,iy), 



dr 



dr^ 



it means that W{r^,Q) = if and only if -^u{r-^) = 0. Consequently our 
result would follow if we can prove that there is no fixed point (^=,,,0) 
= 0. To show this fact, notice first that 



for which -^u^r 



(30) 



d 
dr 



u{r) 



c^(l + Kr^) fi;mr(4Kr^ — 3) 
r2 4(1 - Kr2)3/2 



From the definition of u{r) in (26), the identity u{r^) = is equivalent 
to 

n _ ^^2^l/2 _ mr,{h-AKr1) 
^ *^ Ac\l-Kr^) ■ 

Regarding (1 — Kr'^Y^'^ as (1 — /tr^)^/^(l — Kr^), and substituting the 
above expression of (1 — KrlY^"^ into (30) for r = r,,, we obtain that 



K{AKrl - 3) _^ 1 + /tr^ 
5 — Anr'i 



0, 



which leads to the conclusion that r1 = 5 /2k, < 0. Therefore there is 
no fixed point (r^,, 0) such that -^u^r^ 



the proof. 



0. This conclusion completes 

□ 



9.1. The flow in the (r, z/) plane for k < 0. To study the flow of 
system (25) for /t < 0, we will consider the two cases given by Lemma 



6| namely when system ( 25 ) has no fixed points and when it has exactly 
one fixed point. 
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9.1.1. The case of no fixed points. Since k < 0, and system (25) has 
no fixed points, the function u{r), defined in (26), has no zeroes. But 
lim.r^ou{r) = +00, so u{r) > for all r > 0. Then u{r)/r'^ > for all 
r > 0. Since limr^o g2{r, u) = 0, it follows that linir^o W{r, z/) = +oo, 
where W{r, u) (defined in ([29|) forms the right hand side of the second 
equation in system (25). Since system (25) has no fixed points, W 



IS 



doesn't vanish. Therefore W{r, u) > for all r > and u. 

Notice that the slope of the vector field, u(r, z/), defined in (28) 
of the form v{r,i>) = W{r,i')/i>, which implies that the fiow crosses 
the r axis perpendicularly at every point with r > 0. Also, for r fixed, 
limj,_s.-|-oo W{r, v) = +oo. Moreover, for u fixed, limj._>.oo W{r, v) = 0. 
This means that the fiow has a simple behavior as in Figure ga). These 
orbits correspond to those stated in Theorem ItI (iv). 



r, ' ' '' ' ; ^ . 





— ^ ' ' "■ r i ; ; ' : ] - 'i \ — i ^. \ — r-^ 



V\ ----- 




Figure 4. A sketch of the fiow of system ([25]) for (a) 
K = —2, c = 2, and m = 4, typical for no fixed points; 
(b) K = — 2,c = 2, and m = 6.2, typical for one fixed 
point. 



9.1.2. The case of one fixed point. We start with analyzing the 
behavior of the fiow near the unique fixed point (ro, 0). Let F{r, u) = u 
denote the right hand side in the first equation of system (25). Then 
f F(ro,0) = 0, |;F(ro,0) = 1, and i;W{ro,0) = 0. To determine the 
sign of ^H^(ro,0), notice first that \imr^oW{r, u) = +oo. Since the 
equation W{r, u) = has a single root of the form (ro, 0), with ro > 0, 
it follows that W{r, 0) > for < r < ro. 

To show that W{r, 0) < for r > ro, assume the contrary, which 
(given the fact that ro is the only zero of W{r, 0)) means that W{r, 0) > 
for r > ro. So W^(r, 0) > 0, with equality only for r = ro. But recall 
that we are in the case when the parameters satisfy the inequality 
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c'^K, + > 0. Then a slight variation of the parameters k < 0,m > 0, 
and c 7^ 0, within the region defined by the above inequahty, leads to 
two zeroes for W{r,0), a fact which contradicts Lemma [6] Therefore, 
necessarily, W{r,0) < for r > tq. 

Consequently W{r,0) is decreasing in a small neighborhood of ro, 
so ^iy(ro,0) < 0. But by Lemma [7| ■^W{rQ,0) ^ 0, so necessarily 
W{rQ, 0) < 0. The eigenvalues corresponding to the system obtained 



dr 



by linearizing equations (25) around the fixed point (rg, 0) are given by 
the equation 



(31) 



det 



f^(ro,0) 



1 

-A 



0, 



so these eigenvalues are purely imaginary. In terms of system (25), 
this means that the fixed point (ro, 0) could be a spiral sink, a spiral 
source, or a center. The symmetry of the fiow with respect to the r 
axis excludes the first two possibilities, consequently (ro, 0) is a center 
(see Figure |4](b)). We thus proved that, in a neighborhood of this fixed 
point, there exist infinitely many periodic Eulerian solutions whose 
existence was stated in Theorem [T] (iii) 



To complete the analysis of the fiow of system (25), we will use the 
nuUcline z> = 0, which is given by the equation 



(32) 



1 



c (1 — Kr^ 



m(5 — Ak'^) 
4r2(l - fi;r2)i/2 



Along this curve, which passes through the fixed point (ro,0), and is 
symmetric with respect to the r axis, the vector field has slope zero. 
To understand the qualitative behavior of this curve, notice that 



lim 



K,r 



Kr 



m{5 — 



4r2(l 



Kr 



2)1/2 



But we are restricted to the parameter region given by the inequality 
+ Kc'^ > 0, which is equivalent to 



{-Ky/^ - {-Ky][m{-Ky/^ + {-Ky] > o. 



Since the second factor of this product is positive, it follows that the 
first factor must be positive, therefore the above limit is positive. Con- 



sequently the curve given in (32) is bounded by the horizontal lines 

-.211/2 



z/ = [m(-K)^/2 + KcY''^ and u 



-kV/^ + KC- 



Inside the curve, the vector field has negative slope for u > and pos- 
itive slope for z/ < 0. Outside the curve, the vector field has positive 
slope for z/ > 0, but negative slope for z/ < 0. So the orbits of the 
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flow that stay outside the nullchne curve are unbounded. They corre- 
spond to solutions whose existence was stated in Theorem [T] (iv). This 
conclusion completes the proof of Theorem [7} 

10. Hyperbolic homographic solutions 

In this last section we consider a certain class of homographic orbits, 
which occur only in spaces of negative curvature. In the case K — 1, 
we proved in |2] the existence of hyperbolic Eulerian relative equilibria 
of the curved 3-body problem with equal masses. These orbits behave 
as follows: three bodies of equal masses move along three fixed hyper- 
bolas, each body on one of them; the middle hyperbola, which is a 
geodesic passing through the vertex of the hyperboloid, lies in a plane 
of that is parallel and equidistant from the planes containing the 
other two hyperbolas, none of which is a geodesic. At every moment 
in time, the bodies lie equidistantly from each other on a geodesic hy- 
perbola that rotates hyperbolically. These solutions are the hyperbolic 
counterpart of Eulerian solutions, in the sense that the bodies stay on 
the same geodesic, which rotates hyperbolically, instead of circularly. 
The existence proof we gave in works for any k < 0. We therefore 
provide the following definitions. 

Definition 7. A solution of the curved 3-body problem is called hyper- 
bolic homographic if the bodies maintain a configuration similar to itself 
while rotating hyperbolically. When the bodies remain on the same hy- 
perbolically rotating geodesic, the solution is called hyperbolic Eulerian. 



While there is, so far, no evidence of hyperbolic non- Eulerian homo- 
graphic solutions, we showed in that hyperbolic Eulerian orbits exist 
in the case of equal masses. In the particular case of equal masses, it is 
natural to assume that the middle body moves on a geodesic passing 
through the point (0, 0, (the vertex of the hyperboloid's up- 

per sheet), while the other two bodies are on the same (hyperbolically 
rotating) geodesic, and equidistant from it. 

Consequently we can seek hyperbolic Eulerian solutions of equal 
masses of the form: 

(33) q = (qi, q2, qs), with q^ = {xi, yi, Zi), i = 1, 2, 3, and 

Xi = 0, 7/1 = sinh w, Zi = {kI^^^"^ coshoo, 

^2 = (p^ + /«~^)^^^5 1/2 = psinhcj, Z2 = pcoshcu, 

X3 = — (p^ + K^^)^/^, ?/3 = psinhcj, Z3 = pcoshw. 
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where p := p{t) is the size function and uj := uj{t) is the angular 
function. 

Indeed, for every time t, we have that x1{t) + yf{t) —zf{t) = i = 
1, 2, 3, which means that the bodies stay on the surface H^, while lying 
on the same, possibly (hyperbolically) rotating, geodesic. Therefore 



representation (33 ) of the hypebolic Eulerian homographic orbits agrees 
with Definition [71 

With the help of this analytic representation, we can define Eulerian 
homothetic orbits and hyperbolic Eulerian relative equilibria. 

Definition 8. A hyperbolic Eulerian homographic solution is called 
Eulerian homothetic if the configuration of the system expands or con- 
tracts, but does not rotate hyperbolically. 



In terms of representation (33), an Eulerian homothetic solution oc- 
curs when u{t) is constant, but p{t) is not constant. A straightforward 
computation shows that if uj{t) is constant, the bodies lie initially on 
a geodesic, and the initial velocities are such that the bodies move 
along the geodesic towards or away from a triple collision at the point 
occupied by the fixed body. 

Notice that Definition |8] leads to the same orbits produced by Defini- 
tion |5} While the configuration of the former solution does not rotate 
hyperbolically, and the configuration of the latter solution does not 
rotate elliptically, both fail to rotate while expanding or contracting. 
This is the reason why Definitions |5] and [s] use the same name (Eulerian 
homothetic) for these types of orbits. 

Definition 9. A hyperbolic Eulerian homographic solution is called 
a hyperbolic Eulerian relative equilibrium if the configuration rotates 
hyperbolically while its size remains constant. 



In terms of representation (33), hyperbolic Eulerian relative equilib- 
ria occur when p{t) is constant, while uj{t) is not constant. 

Unlike for Lagrangian and Eulerian solutions, hyperbolic Eulerian 
homographic orbits exist only in the form of homothetic solutions or 
relative equilibria. As we will further prove, any composition between 
a homothetic orbit and a relative equilibrium fails to be a solution of 
system 



Theorem 8. In the curved 3-body problem of equal masses with k, < 0, 
the only hyperbolic Eulerian homographic solutions are either Eulerian 
homothetic orbits or hyperbolic Eulerian relative equilibria. 



32 



Florin Diacu and Ernesto Perez-Chavela 



Proof. Consider for system ([8j) a solution of the form ( 33 ) that is not 
homothetic. Then 

1/2, 



/«qi q2 = /«qi qs = ' p, 



/«q2 qs = -1 - 2kp 



Xi 



xi = 0, 



yi 



\K\ 



'^I'^Oj cosh 



'^^"^ sinhw, 



X2 = -xs 



PP 



Kq2 q2 



Z/2 
.2 



(p2 + fi:-i)i/2' 
Vs = p sinh u + pu cosh w, 
^3 = pcoshw + pwsinho;, 



Kq2 q2 = /tqs qs 



2 • 2 



Kp 



.2 ' 



a^2 = -a^s 



(p2 + K-l)l/2 (p2 + ^-l)S/2' 

i/2 = Z/s = (p + pcj^) sinh w + (pw + 2pLo) cosh w, 
^2 = ^s = (p + p(^'^) cosh w + (pci) + 2puj) sinh cu. 
Substituting these expressions in system (|8]), we are led to an identity 
corresponding to Xi. The other equations lead to the system 



X2,X3 : 


E = 




yi ■ 


K ^-^/^a) cosh a; = 0, 




zi : 


\ K,\~^^'^ Co sinh u = 0, 




1/2, i/s : 


E sinh + F cosh 


= 0, 


^2, ^S : 


E cosh + F sinh a; 


= 0, 



where 



F:=F(t)=p + p(l + .pV-T^+ "^(^-^'^^^ 



1 + Kp2 4p2|l + Kp2|l/2^ 

F := F(t) = pu + 2puj. 
This system can be obviously satisfied only if 

'u = 

(34) <; a) = - 



2pcj 
P 



The first equation implies that uj{t) = at + b, where a and b are con- 
stants, which means that co(t) = a. Since we assumed that the solution 
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is not homothetic, we necessarily have a 7^ 0. But from the second 
equation, we can conclude that 



where c is a constant. Since a 7^ 0, it follows that p(t) is constant, 
which means that the homographic solution is a relative equilibrium. 
This conclusion is also satisfied by the third equation, which reduces 
to 



9 r, mil — 4:Kp) 



4p3|l + Kp2|l/2' 

being verified by two values of a (equal in absolute value, but of oppo- 
site signs) for every k and p fixed. Therefore every hyperbolic Eulerian 
homographic solution that is not Eulerian homothetic is a hyperbolic 
Eulerian relative equilibrium. This conclusion completes the proof. □ 

Since a slight perturbation of hyperbolic Eulerian relative equilibria, 
within the set of hyperbolic Eulerian homographic solutions, produces 
no orbits with variable size, it means that hyperbolic Eulerian rela- 
tive equilibria of equal masses are unstable. So though they exist in a 
mathematical sense, as proved above (as well as in [2], using a direct 
method), such equal- mass orbits are unlikely to be found in a (hypo- 
thetical) hyperbolic physical universe. 
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